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The power spectrum of quantum fluctuations of the electromagnetic field produced by an elemen- 
tary particle is determined. It is found that in a wide range of practically important frequencies the 
power spectrum of fluctuations exhibits an inverse frequency dependence. The magnitude of fluctu- 
l/") ' ations produced by a conducting sample is shown to have a Gaussian distribution around its mean 

value, and its dependence on the sample geometry is determined. In particular, it is demonstrated 
that for geometrically similar samples the power spectrum is inversely proportional to the sample 
volume. It is argued also that the magnitude of fluctuations induced by external electric field is 
proportional to the field strength squared. A comparison with experimental data on flicker noise 
measurements in continuous metal films is made. 
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I. INTRODUCTION 



Measurements of voltage fluctuations in various media show that at sufficiently low frequencies, power spectra of 
fluctuations in all conducting materials exhibit a universal profile which is close to inverse frequency dependence, and 
called for this reason a flicker l//-noise. Although this noise is dominating only at low frequencies, experiments show 
the presence of the l//-component in the whole measured band, from 10 -6 Hz to 10 6 Hz. Also, the following three 
main characteristic properties are universally held: The l//-noise produced by a conducting sample 1) is inversely 
O proportional to its volume, 2) is Gaussian, and 3) its part induced by external electric field is proportional to the field 
l— "™ ', strength squared. 

' A number of mechanisms has been put forward to explain the origin of the l//-noise. The property 3) is typical 



for resistive systems, and suggests that the flicker noise can arise from resistance fluctuations pj. It has been also 
I/"") ■ proposed that fluctuations in the carrier mobility 0,0, or in the number of carriers caused by surface traps Q might 
explain the origin of the l//-profile of the power spectrum. All these models, however, have restricted validity, because 
they involve one or another assumption specific to the problem under consideration. For instance, assuming that the 
resistance fluctuations in the first of the above-mentioned models spring from temperature fluctuations, one has to 
choose an appropriate spatial correlation of these fluctuations in order to obtain the l//-profile of the power spectrum 
Similarly, the model suggested in Ref. £| requires a specific distribution of trapping times, etc. In addition to that, 
the models suggested so far reproduce the l//-profile only in a restricted range of frequencies. On the other hand, the 
ubiquity of the flicker noise, and universality of its properties imply that there must exist a simple and universal, and 
therefore, fundamental underlying reason. It is natural to look for this reason in the quantum properties of charge 
carriers. In this direction, the problem has been extensively investigated by Handel and co-workers [(|. Handel's 
approach is based on the theory of infrared radiative corrections in quantum electrodynamics. Handel showed that 
O the 1// power spectrum of photons emitted in any scattering process can be derived from the well- know property 
of bremsstrahlung, namely, from the infrared divergence of the cross-section considered as a function of the energy 
• i-h , loss. Thus, this theory treats the l//-noise as a relativistic effect (in fact, the noise level in this theory ~ a(Av) 2 J c 2 , 
^ ■ where a is the fine structure constant, Av velocity change of the particle being scattered, and c the speed of light). 
It should be mentioned, however, that the Handel's theory has been severely criticized in many respects 0,0- 

The purpose of this paper is to draw attention to another quantum aspect of the electromagnetic interaction, 
which is of a purely nonrelativistic nature. It turns out that there is a simple and quite general property of the 
electromagnetic interactions of quantized matter, which may be the origin of flicker noise. Namely, it will be shown 
below that the power spectrum of the Coulomb field fluctuations produced by spreading wave packet of a free charged 
particle exhibits an inverse frequency dependence in the low- frequency limit. We will show also that the power 
spectrum possesses the three above-mentioned characteristic properties of flicker noise as well. 
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The main tool we use in investigating quantum fluctuations of the Coulomb field produced by a quantum particle 
is the two-point correlation function of the electromagnetic field. Properties of this function in the coincidence limit 
were investigated in detail in Ref. It was found, in particular, that the root mean square fluctuation of the 
Coulomb potential is of zeroth order in the Planck constant h. It was shown, furthermore, that the Fourier transform 
of the correlation function with respect to time exhibits an inverse frequency dependence in the low- frequency limit. 
One of the goals of the present paper is to show that essentially the same formula describes also the low-frequency 
limit of the power spectrum of fluctuating electromagnetic fields measured at two distinct time instants. The reason 
underlying this similarity is that in both cases the leading contribution to the correlation function is contained in its 
disconnected part. 

The paper is organized as follows. In Sec.[n]a preliminary consideration of the problem is given, and the Schwinger- 
Keldysh formalism used throughout the work is briefly reviewed. The low-frequency asymptotic of the power spectrum 
is calculated in Sec. IIIII It is proved in Sec. IIII Al that the low-frequency asymptotic of the connected part of the 
correlation function is logarithmic. Contribution of the disconnected part is calculated exactly in Sec. IIII Bl and is 
shown to exhibit an inverse frequency dependence. The obtained result is analyzed to verify that it is in agreement with 
the above-mentioned general characteristics of flicker noise, and then compared with experimental data in Sec. IIII CI 
Section llVl summarizes the results of the work and states the conclusion. 



II. PRELIMINARIES 



Let us consider a single particle with mass m and electric charge e. In classical theory, the electromagnetic field 
produced by such particle at rest is described by the Coulomb potential 

^° = i^ a = o. a) 

Our aim is to determine quantum properties of this potential. More precisely, we will be concerned with the quantum 
fluctuations of the field-theoretic counterpart of this potential at the zeroth order in the Planck constant h. 

In quantum theory, Eq. is reproduced by calculating the corresponding mean fields, (in|Ao|in), (in|A|in). In the 
Schwinger-Keldysh formalism fLUL UA\ , they are given by 

(mli^in) = J D$_ J D$+ A^exp{iS[$+\ -iS[$-]}, (2) 

where the subscript + (— ) shows that the time argument of the integration variable runs from — oo to +oo (from +oo 
to — oo). Integration is over all fields satisfying 

4>t^<t>t, for t -> -oo, 

$ + = for t -> +oo , (3) 

where $ collectively denotes the fundamental fields of the theory, $ = </>*}, the components 4>,(f)* describing 

the charged particle which for simplicity will be assumed scalar, (j)Q is the given particle state, D<S> the invariant 
integral measure, and S — S[<f>] the action functional of the system. Assuming that the gradient invariance of the 
theory is fixed by the Lorentz condition 

G = 0M„ = O, (4) 

the action takes the form 

5[*] = S ffl + Sgf[$], 

S [$] - Jd 4 x{(d^* +ieA^*){d^~ieA^)-m 2 ( f ) *^}-^Jd 4 xF^F^ , 
Sgf[9] = -\J dAx G 2 > F ^ = d » A » - dvAp . (5) 

Diagrammatics generated upon expanding the integral 121 in powers of the coupling constant e consists of the following 
elements. There are four types of pairings for each field A^ or 0, corresponding to the four different ways of placing 
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two field operators on the two branches of the time path. They are conveniently combined into 2x2 matrices 1 



W(x,y) 



D{x,y) 



Dl u f (x,y) D?_(x,y) \ = i{TA»{x)A»(y)) i(A»(y)A»(x)) 
Dt + (x,y) DfL{x,y) J { i(A»(x)A»(y)) i(fA»(x)A»(y)) 



D ++ (x,y) D + -(x,y)\ = (i{T$(x)ft(y)) i{ft(y)$(x)} 
D-+{x,y) D — (x,y) J \ i(fa)ft(y)) i(f$(x)4>Hy)) 



where the operation of time ordering T (T) arranges the factors so that the time arguments decrease (increase) from 
left to right, and (-)o denotes vacuum averaging. The "propagators" D satisfy the following matrix equations 

d A z®^(x,z)V™(z,y) = -cS^(x-y), <V (s, y) = i g^Ig^ry (6) 

rf 4 z <S(x,z)Q(z t y) = -cS^(x-y), <S(x,y)=i J! S f , (7) 

d(j)*(x)d(j){y) 

where is the free field part of the action, e, i are 2x2 matrices with respect to indices +, — : 

1 0\ . _ f 1 
1 J ' l ~ ^ -1 

As in the ordinary Feynman diagrammatics of the S-matrix theory, the propagators are contracted with the vertex 
factors generated by the interaction part of the action, 5 mt [$] = — S , ' 2 - ) [$], with subsequent summation over 
(+,— ) in the vertices, each "— " vertex coming with an extra factor (—1). This can be represented as the matrix 
multiplication of 2)p„,£) with suitable matrix vertices. For instance, the Adcjxjf part of the action generates the 
matrix vertex 03 which in components has the form 

V iik( x >V^ z ) = s vk 



where the indices k take the values +, — , and Sijk is defined by s+++ = —s = 1 and zero otherwise. External 

(f> (</)*) line is represented in this notation by a column (row) 



satisfying 

J d A z <8(x,z)t(z) = ( , J d A x r + «) = (0,0). (8) 

The tree diagrams contributing to the right hand side of Eq. 10 are depicted in Fig. ^ Of these only the diagram 
in Fig. ^a) gives rise to a non-zero contribution. The momentum flow in the other diagram is inconsistent with the 
momentum conservation in the vertex for any non-zero value of the momentum transfer p^, because the lines of this 
diagram are all on the mass shell. The same diagram^a) represents also the tree value of the in-out matrix element 
(out|A A '(a;)|in) calculated using the ordinary Feynman rules. This is natural because we are concerned here with 
a one-point function whose matrix element is evaluated between one-particle states for which |in) is essentially the 
same as |out). The things change, however, if we want to determine the field fluctuation. In the Schwinger-Keldysh 
formalism, the expectation value of the product A^{x)A u {x') has the form 

{m\A»(x)A v (x')\m) = J D$_ J D$+ A»(x)A u + {x')exp{iS[$ + ]-iS[®-]}, (9) 



1 Below, Gothic letters are used to distinguish quantities representing columns, matrices etc. with respect to indices 
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The corresponding tree diagrams are shown in Fig. [5] We see that these diagrams differ from those we would have 
obtained applying the ordinary Feynman rules to the quantity (out|A M (x)A v (x')|in) . We define the correlation function 
of the electromagnetic field fluctuations by 

C^(x,x') = (m\A^x)A u (x')\m) - (m\A^x)\m)(m\A u (x')\m). (10) 

The first term in this definition is usually written in a symmetrical form thus rendering the function C tiv {x,x') real. 
However, this makes no difference for what follows because it will be shown below that the leading low-frequency term 
of the correlation function is contained entirely in its disconnected part [the second term in Eq. p()|l]. 

We are concerned with the power spectrum of correlations in the values of the electromagnetic fields measured at 
two distinct time instants (spatial separation between the observation points, \x — x'\, is also kept arbitrary). Thus, 
fixing one of the time arguments, say, t', we define the Fourier transform of C^ v {x^x r ) with respect to (t — t') 

C flt/ {x,x',t',u;) = { dtC^&x'^-t'K (11) 



To complete the present section, we write out explicit expressions for various pairings of the photon and scalar fields 

f d 4 k e - ik ( x -v) f d 4 k e -ik(x-y) 



(2tt) 4 to 2 - k 2 - iO ' ~ v ' J (2tt) 4 k 2 - m 2 - iO ' 
D-+U-U) = > I 7^3 6(k°)6(k 2 - m 2 )e- lk(x -y*> , D + -{x, y) = D_+(y, x) . (12) 



III. EVALUATION OF THE LEADING CONTRIBUTION 



Evaluation of the low-frequency asymptotic of the correlation function proceeds in two steps. First, we will prove 
in Sec. IIII Al that the low-frequency asymptotic of the connected part of C M „ [the first term in Eq. IjlUfl ] is only 
logarithmic. The contribution of the disconnected part will be calculated exactly in Sec. IIII 61 It will be shown that 
this contribution exhibits an inverse frequency dependence, and thus dominates in the low- frequency limit. 



A. Low-frequency asymptotic of the connected part of correlation function 

Before going into detailed calculations, let us first exclude the diagrams in Fig. [5] which do not contain the h° 
contribution. It is not difficult to see that such are the diagrams without internal matter line. Indeed, consider, for 
instance, the diagram (g). It is proportional to the integral 

/ p ik{x~x ) 
d 4^ (fc )(5(fc2) ___ 

which does not involve the particle mass at all. Taking into account that each external scalar line gives rise to the 
factor (2e q ) _1 / 2 , where e q = \J m 2 + q 2 w m , we see that the contribution of the diagraming) is proportional to 
X/m. Hence, on dimensional grounds, this diagram is proportional to h. The same is true of all other diagrams without 
internal matter lines. 

The sum of diagrams (a)-(e) in Fig. [21 has the following symbolic form 

V + 1% , V = 7 [t f 2J a !D3^t] , 

where the superscript "tr" means transposition of the indices and spacetime arguments referring to the points of 
observation: fi <-+ v, + , x <-» x' [the transposed contribution is represented by the diagrams collected in part (e) 
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of Fig. [5] ■ Written longhand, reads 



I liv {x,x')=ie 2 II d 4 zdVj + D° ++ (x,z) 

- D° ++ (x,z) 

- D° + _(x,z) 

+ D°,(x,z) 



<ft(z)d li D ++ (z,z r )dlMz') 



D° _(*',*') 
D°__(z',x') 
D° + _(z',x') 

D°__(z',x')\, 



where 



(13) 



Contribution of the third term in the right hand side of Eq. <|13[1 is zero identically. Indeed, using Eq. (|12fl . and 
performing spacetime integrations we see that the three lines coming, say, into z-vertex are all on the mass shell, 
which is inconsistent with the momentum conservation in the vertex. Taking the Fourier transform of I^ u {x, x'), 



I„ v (x,x',t , ,u) = / dtl^x.x'Y^t-t') 



the remaining terms in Eq. (|13|) take the form 



d 3 q d 3 p a*(q)a(q + p 



(27r)3(27r)3 



2 e -i P °{t'-t )+i P x' J (pj gx _ aj / ) w j 



(14) 



q+p 



-q+p tq ; 



where 



J^ip.q.x- x',oj) = -i 



d 3 k 

(2^ 



x {D a ++ (k)D ++ (q + k)D\_(k -p)- D\ + {k)D + _{q + k)D°_ _(k - p) 
+D° + _(k)D^(q + k)D°__(k - p)} kQ ^ 



(15) 



Here e q = \/ q 2 + m 2 is the particle's energy, to, e, g M , and a(q) are its mass, electric charge, 4-momentum, and 
momentum wave function at some time instant to, p M 4-momentum transfer. The function a(q) is normalized by 



d 3 q 



>(q)| 2 = l 



and is generally of the form 



a(q) = b(q)e 



tqx 



(16) 



(17) 



where x is the particle mean position, and b(q) describes the momentum space profile of the particle wave packet. 
For simplicity, the particle charge distribution will be assumed spherically-symmetric in what follows. Then b(q) is a 
function of q 2 : b(q) = f3(q 2 )- 

As in Ref. @, we work within the long-range expansion which implies, in particular, that the spatial separations 
between the field-producing particle (more precisely, its mean position) and the points of observation, r = \x — xq\, 
r' = \x' — xq\, are such that rq ^> 1, r'q ^> 1, where q — -J (q 2 ) is the particle momentum variance (note that q is 
time-independent for a free particle). To the leading order of the long-range expansion, p^, k^ in the vertex factors 
in Eq. (|15f) can be neglected in comparison with q^. 
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Let us now show that the low- frequency asymptotic of is weaker than 1/ui. The first and the third terms in the 
integrand in Eq. I|15f) give rise to a contribution which is finite at uj = 0. We have for the third term 



J d 3 k e tk( > x - x ^D a + _(k)D—(q + k)D __(k-p) 



k°=0 



d 2 o / d\k 



HI 



[-2(qk) - iO] [p 2 + 2{pk) - i0] A\q\p 2 J cos 9 + iO 



2\q\p 2 



Similarly, 



J d 3 k e ik ^-^D° ++ (k)D ++ (q + k)D° + _(k - p) 



fc°=0 



d 2 o J d\k 

o 



k 2 e lk ( x - x ">2iri9{p°)5 ((p ) 2 -(k-p) 
[fc 2 - i0}[k 2 + 2(qk) - iO] 



is also finite for all p and q ^ [the singularity at q = or p = in these expressions is inessential as it is removed 
by the factors q 2 and p 2 in the integral measure in Eq. I|14|l ]. The second term in the integrand, however, leads to a 
divergence at uj = 0. This divergence comes from integration over small |fc|. Hence, it is the same as the divergence of 
the integral 



dfo j d\k 

o 



k 2 5{uj 2 + 2mu -k 2 2{qk)) 
[-oj 2 + k 2 - iO][p 2 - iO] 



y/q 2 +u 2 +2mu+\q\ 



\q\p 2 



udu 



iO 



■sj q 2 +u> 2 +2mLO—\q\ 



The latter integral diverges for uj — > only logarithmically. 



B. Low- frequency asymptotic of the disconnected part of correlation function 

Let us turn to the disconnected part of the correlation function. To find its Fourier transform, we have to evaluate 
the following integral 

*(in|A M (x)|in)e fa '(*-*'> = q„e- iut ' I{r,w) , (18) 



where 



Substituting 



gives 



-t-oo 



r = x — Xq . 



2m 2m 



For a spherically-symmetric charge distribution, b(q) = f3(q 2 ), I is actually a function of r = \r\, I(r,ui) — I(r,u>), 
and hence, averaging over directions of r, one can write 



to f d 3 q 2 2 „ ? sin 6> sulfur) 



27rr 7 (2tt) 3 ^ u^q 2 cos 2 g + 2muj 



- \q\ cos 9 + \/q 2 cos 2 6> + 2muj , (19) 
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where 9 is the angle between the vectors p, q , and it is assumed that to > 0. Taking u as the integration variable 
yields 

+oo V q 2 +2mw+|q 

me M *° /" 9 . , r,. . 9 , /" , sm(itr) , „ 

= "g^T / dq 2 p*(q 2 )(3(q 2 + 2mw) I du ^T^ ■ ( 20 ) 

% /q 2 +2mw-|q| 

To further transform this integral, it is convenient to define a function T(q 2 ,uj) according to 

T(q 2 ,uj) = ^ J dzf3*(z)f3(z + 2mu). (21) 

q 2 

Then integrating by parts, and taking into account that r(<7 2 ,o;) — > for q 2 — > oo, brings Eq. (|2(J|) to the form 

+ 00 

~ me^Q f r(q 2 ,^) f sin(w) 



+ 



2?rr J V? 2 + 2mw 
o 

sin(itr) 

u—^J q 2 +2muj—\q\ 



l—^J q 2 +2maj+|q| 

(22) 



Finally, integrating / by parts once more, we find 

+ oo 

J(r,w) = y d|«| sin (V q 2 + 2muj ^2T cos(|q|r) + -— sin(|q|r)| . (23) 

o 

This exact expression considerably simplifies in the practically important case of low to and large r. Namely, if u is 
such that 

lo < — = uj , (24) 
mr 

and also 

r<? > 1 

(and therefore, oj <C q 2 /m), then the first term in / turns out to be exponentially small (~ e~ r? ) because of the 
oscillating product of trigonometric functions. Replacing sin 2 (gr) by its average value (1/2) in the rest of / gives 

~\~ OO ~\~ OO 

I(r,u) = / d\q\-^- = - e - ' j = ' 2 / dz\p{z)\ 2 . (25) 

47rr 2 cj J d\q\ Anr 2 uj 167r 3 r 2 u; J 1 v ;l v y 

o o 

Substituting the obtained expression into the defining equations (|1H|) . I|10|l . JUJ, we thus obtain the following expres- 
sion for the low-frequency asymptotic of the correlation function 

2 3 3 ( ' ) 



(26) 

all other components of C MJ , being suppressed by the factor \q\/m <C 1. For the time instants if such that r, r' ^> 
where D t i is a characteristic linear dimension of the particle charge distribution (for instance, variance of the particle 
coordinates), Eq. (12(it simplifies to 

C 00 (x,x',t',uj) = e ^(to-t') f f dz \^ z )\\ (27) 
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In applications to microelectronics, ui varies from 10~ 6 Hz to 10 6 Hz, the relevant distances r are usually 10 -5 cm to 
10~ 2 cm, q ~ h/d, where d « 10~ 8 cm is the lattice spacing, and m is the effective electron mass, m « 10 _27 g, hence, 
ujq w I0 10 Hz, so the conditions r 3> d, u> <C wo are always well-satisfied. 

Equation l|27l) represents an individual contribution of a conduction electron to the electric potential fluctuation. 
Considering a large number of uncorrelated electrons in a sample, one should take into account that the corresponding 
time instants to are distributed uniformly. Because of the oscillating exponent e 41 ^* 0- * ), the magnitude of the total 
noise remains at the level of the individual contribution independently of the number of electrons. Therefore, summing 
up all contributions amounts simply to averaging over Xq : 



\C$(x,x',u;)\ = — ^— / dz\f3{z)\ 2 I ^ , (28) 



+ 00 

QAtT A \u\VL 

where f2 is the sample volume. To visualize it, one can say that various uncompensated individual contributions 
"flicker" in various points of the sample. 

It is clear from the above discussion that the distribution of the noise magnitude around the value Cqq* is Gaussian, 
by virtue of the central limiting theorem. Thus, the quantum field fluctuations in a sample possess the above-mentioned 
property 2). 

Next, let us consider the case when the sample is in an external electric field, E. Under the influence of the external 
field, both the electron wave functions and the statistical probability distribution are changed. It turns out, however, 
that in practice, the change of the wave functions is negligible. Namely, a direct calculation shows that the relative 
value of the correction is of the order 

h 2 me\E\ 



q±L ' 

where L is the smallest linear dimension of the sample. Using the above estimates for q, L in the condition n q -C I 
gives \E\ <C 10 10 , in the cgs system of units. In applications to microelectronics, the field strength is usually \E\ < 1, 
so K q « 10~ 10 . On the other hand, the relative change of the statistical distribution function, given by the kinetic 
theory, is of the order 

e\E\l 



where I is the electron mean free path, and e its mean energy, e w ti 2 /md 2 . n s is generally not small. Furthermore, 
the statistical distribution is generally analytic in E in a vicinity of E = 0, and hence the leading correction to the 
scalar quantity Coo is proportional to E 2 [property 3)]. 

Because of collisions with phonons and impurities in the crystal, evolution of the conduction electrons usually cannot 
be considered free. It is important, however, that in view of the smallness of the electron mass in comparison with 
the mass of the lattice atoms, these collisions may often be considered elastic. In such collisions, the electron energy 
is not changed, and therefore, so is time evolution of the electron wave function. Hence, the above results concerning 
dispersion of the electromagnetic field fluctuations remain essentially the same upon account of the electron collisions, 
except that now they must be expressed in terms of the electron density matrix, rather than wave function. Denoting 
the momentum space density matrix of the electron by p(p 2 , q 2 ), Eq. I|28|) thus takes the form 

+00 

n 

The quantity denoted here by To is nothing but the expectation value of the inverse particle momentum: 

, d 3 q p(q 2 ,q 2 ) / 1 

1 — 



(2tt)3 \q\ \\q 

Hence, by the order of magnitude, r ~ 1/q (in the ordinary units, T ~ h/q.) 

In practice, one is interested in fluctuations of the voltage, U, between two leads attached to a sample. Using the 
above results, it is not difficult to write down an expression for the voltage correlation function, C\j(x,x',cj). We 
have 

(in\U(x,x',t)\in)(m\U(x,x',t')\m) = (m\A Q (x,t) - A (x',t)\in)(in\Ao(x,t') - A (x\t')\m) 
= (in|A (a;,i)|in)(m|Ao(a;,< , )|in) + (in\A (x' ,t)\m)(m\A (x' ,t')\in) 

-(in|^ (x / ,t)|in)(in|ilo(aJ,ON - HA Q (x,t)\m)(iii\A (x' ,t')\in) . (30) 
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Fourier transforming and applying Eq. 127|) to each of the four terms in this expression, we obtain the low-frequency 
asymptotic of the overall voltage fluctuation across the sample 



Ctot 
i U 



where 



G = 



(x, x 1 , 



e 2 r G 
16ir 2 \u\n 



,...11 1 1 



(31) 



(32) 



Let us discuss the role of the sample geometry in somewhat more detail. Consider two geometrically similar samples, 
and let n be the ratio of their linear dimensions (see Fig- - Such samples are characterized by the same value of 
the G-factor. Indeed, the voltage across each sample is measured via two leads attached to its surface. Hence, the 
radius-vectors of the leads drawn from the center of similitude scale by the same factor n, and therefore, 



d 3 



d 3 x 



n 3 ni 



\nx — Xo\ 2 \nx' — Xq 
cPxq 



Xq\ 



Wixo) 



\nx — nxo\ 2 \nx' — uxq\ 



d 3 x 



and likewise for the other terms in Eq. I|32|) . Thus, it follows from Eq. (|28[1 that for geometrically similar samples, the 
noise level is inversely proportional to the sample volume [property 1)]. 

It is worth also to make the following comment concerning Eq. (|32Jl . The integrand in this formula involves the 
terms 1 jr 3 and 1 jr' 3 which give rise formally to a logarithmic divergence when the observation points approach the 
sample. In this connection, it should be recalled that the above calculations have been carried out under the condition 



rg > 1, hence, 



cannot be taken too small. Furthermore, one should remember that in any field measurement 



in a given point, one deals actually with the field averaged over a small but finite domain surrounding this point, 
i.e., the voltage lead in our case. Thus, for instance, the quantity Ao(x,t) appearing in the expression l|3U|l is to be 
substituted by 



where Dl is the voltage lead volume. As the result of this substitution, the term 1/r 3 takes the form 

1 f f d 3 xd 3 x 



n 2 



x \ 2 \x 



This is finite however the lead is attached to the sample. Equation l|32|l is recovered in the case when r, r' are much 
larger than the sample size. 



C. A preliminary comparison with experiment 

For a preliminary comparison of the noise level given by Eq. I|31f) with experimental data, let us use the results of the 
classic paper by Voss and Clarke ;5] where the l//-noise in continuous Bi and Au films was measured. The bismuth 
and gold samples had dimensions lOOOAxlO^m x 120 /im and 250Ax8/im x 625 /im, respectively. The spectrum 
power for / = lo/2tt = 1Hz and \E\ — was found to be about 10~ 14 V 2 / H z in the case of Bi sample, and about 
10~ 16 V 2 / Hz in the case of Au sample (it is assumed here that the contributions from external pickup and amplifier 
noise in the experimental setup of Ref. [5| do not dominate over the sample noise component). 

For a rough estimate of the right hand side of Eq. (|31|l we take Tq w H/q ~ d w 10 -8 cm, and note that fl w 10 -10 cm 3 
for both samples. However, the G-factor takes different values for the two sample configurations. It is seen from 
Eq. (|32|l that for elongated samples of fixed volume, G decreases roughly as the inverse sample length. Evaluating the 
G-factor as explained in the preceding section gives G = 0.03 for Bi sample, and G = 0.007 for Au sample. Substitution 
of the above estimates into Eq. i|3~T|) gives the order of magnitude of the noise spectrum power 10~ 15 V 2 /Hz, and 
10 _16 V 2 /Hz for Bi and Au sample, respectively. We see that in both cases the theoretical estimates are in a 
reasonable agreement with the experimental data. A more detailed comparison is a much more difficult task which 
would require specification of the precise form of the function (3{z). 
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IV. CONCLUSIONS 



We have shown that in the whole measured frequency band, the power spectrum of quantum fluctuations of the 
Coulomb field exhibits an inverse frequency dependence [Cf. Eq. i|od|) ]. This result is obtained by evaluating the two- 
point correlation function of the electromagnetic field within the long-range expansion. It was proved, in particular, 
that the low- frequency asymptotic of the Fourier transform of this function is dominated by its disconnected part. 
Furthermore, we have argued that the derived power spectrum possesses the other experimentally observed properties 
of the 1//- noise as well. In particular, the exact dependence of flicker noise on the sample geometry has been 
established. That the power spectrum of flicker noise is proportional to the sample volume only very approximately 
was noticed long ago, but precise form of this dependence has not been inferred. We have shown that for geometrically 
similar samples the noise power spectrum is inversely proportional to the sample volume, while for samples of the 
same volume dependence on the sample geometry is described by the dimensionless G- factor given by Eq. (|32|l . 

Finally, according to the semi-quantitative estimates of Sec. lIII d the found noise level is in a reasonable agreement 
with experimental data. Thus, quantum fluctuations of the Coulomb field produced by elementary particles can be 
considered as one of the underlying mechanisms of the observed 1 / /-noise. 
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Figure captions 

Fig.l: Tree contribution to the right hand side of Eq. Wavy lines represent photon propagators, solid lines the 
scalar particle. 

Fig. 2: Tree contribution to the right hand side of Eq. @. Part (e) of the figure represents the "transposition" of 
diagrams (a)-(d) (see Sec. IIII B|) . 

Fig. 3: Voltage measurement in geometrically similar samples (n = 2). 
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